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We report on a search for the charmless decays B + — > <f>-K + and B° — > <f>-K° that are strongly 
suppressed in the Standard Model. The analysis is based on a data sample of 657 x 10 B BB pairs 
collected at the T(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. We find no significant signal and set upper limits of 3.3 x 10" 7 for B+ -> 0tt+ and 
1.5 x 1CT 7 for B° -t 4>-k° at the 90% confidence level. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



In the Standard Model (SM), the charmless two-body 
hadronic decays B + —} <pTr + [l[ and B a — > 07r° are 
highly suppressed since they are forbidden at tree level 
and are only possible through the penguin process shown 
in Fig. QJa). The expected SM branching fractions for 
these decays are B(B + -» <fiir + ) ~ 3.2 x 10~ 8 and 
B(B° -> (j)TT ) ~ 6.8 x 10" 9 0, in which the largest 
contribution comes from radiative corrections and uj-ip 
mixing. In some New Physics (NP) scenarios such as 
models with a Z' boson p|, H[ or the Constrained Mini- 
mal Supersymmetric Standard Model (CMSSM) |[, the 
branching fractions could be enhanced up to the 1CP 7 
level. Figure QJb) shows a typical CMSSM contribution 

tO B — > (f)TT. 

Since B — > 07r decays are very sensitive to NP, mea- 
surements of these decays may constrain and potentially 
reveal such contributions. Furthermore, measurements of 
B — » 07T decays also provide a means to study SM con- 
tributions from suppressed diagrams in other important 
decay modes such as B° —> (j)K° [||. A previous search 
by the BaBar collaboration set upper limits of B(B + — > 
07T+) < 2.4 x 1CT 7 and B(B° -+ 0tt o ) < 2.8 x 1CT 7 
at the 90% confidence level (CL) A later measure- 
ment of the three-body inclusive branching fraction for 
B + — > K + K~n + [1] also did not report any evidence for 

B+ -> 07T+. 

In this paper, we report on a search for B + — > cf)ir + 
and B a — > 4>ir° based on a 605 fb _1 data sample, which 



s s 




FIG. 1: (a) The SM three-gluon hairpin penguin diagram 
for B — > 4>tt decays, (b) One of the CMSSM diagrams that 
contributes to B — > <j)n. In both a) and b), the ss quark pair 
hadronizes as a 4> meson. 



corresponds to (657±9) x 10 6 BB events. The data were 
collected with the Belle detector [9] at the KEKB [Io[ 
asymmetric-energy e + e~ collider operating at the T(4S I ) 
resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. 
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An iron flux-return located outside of the coil is instru- 
mented to detect mesons and to identify muons. 
Two inner detector configurations were used: a 2.0 cm 
beampipe and a 3-layer silicon vertex detector were used 
for the first sample of 152 x 10 6 BB pairs, while a 1.5 
cm beampipe, a 4-layer silicon detector, and a small-cell 
inner drift chamber were used to record the remaining 
505 x 10 6 BB pairs [HI- 

To search for B + — > 4>ir + and B° — > 4>ir°, we combine 
<f> — > K + K~ candidates with either a tt + or ir° — > 77. 
Particle identification (PID) for charged kaons from the 
(j> decays and the charged pion is based on the likelihood 
ratios Rk,-k = Lk+l ' wnere Lk and L n denote, respec- 
tively, the individuaf likelihoods for kaons and pions de- 
rived from ACC and TOF information and dE/dx mea- 
surements in the CDC. The PID selections, Rk,tt > 0.3 
for kaon candidates and Rk,k < 0.2 for pion candi- 
dates, are applied to all charged particles. The PID 
efficiencies are 87% (86%) for kaon pairs (high momen- 
tum single pions) in B + (f>n + and 86% for kaon pairs 



in B° 



while the probability of misidcntifying 



a kaon as a pion (a pion as a kaon) is 6% (12%) for 
both modes. Candidate 7r°'s are reconstructed from 7 
pairs that have invariant mass between 115.3 MeV/c 2 
and 152.8 MeV/c 2 , corresponding to ±2.5cr standard de- 
viations (ct). In addition, these photons are required 
to have energies greater than 0.2 GeV. A K + K~ pair 
is required to have an invariant mass within the range 
1.008 GeV/c 2 < M K+K - < 1.031 GeV/c 2 (±2.5 times 
the <fi full width). 

B meson candidates are identified with two kine- 
matic variables: beam-energy-constrained mass, M bc = 
V^bcam - IEiPil 2 > and energy difference AE = 
]C,- -Ej — Sboam, where E beam is the beam energy, and 
pi and Ei are the momenta and energies, respectively, of 
the daughters of the reconstructed B meson candidate 
in the center-of-mass (CM) frame. We fit B candidates 
that lie within the fit region defined by \AE\ < 0.1 GeV 
and M bc > 5.20 GeV/c 2 for B+ -> 0tt+ and \AE\ < 0.4 
GeV and M hc > 5.20 GeV/c 2 for B° -> (jm° . The sig- 
nal regions are defined by \AE\ < 0.04 GeV (±3. Oct) and 
M hc > 5.27 GeV/c 2 (±3.0ct) for B+ -> 0tt+, and -0.16 
GeV (5.0ct) < AE < 0.10 GeV (3.0ct) and M hc > 5.27 
GeV/c 2 (±3. Oct) for B° — > 07r°. We select an asymmet- 
ric signal region for B — > tfiir since photons may inter- 
act with the intervening detector material before entering 
the ECL and there may be energy leakage from the ECL 
crystals. 

The main background arises from the continuum pro- 
cess, e + e~ — > qq, where q = u,d, s, c. To suppress this, 
observables based on the event topology are utilized. The 
event shape in the CM frame is spherical for BB events 
and jet-like for continuum events. This difference is ex- 
ploited by the event-shape variable, which is a Fisher dis- 
criminant formed out of 16 modified Super Fox- Wolfram 
moments [IH, [l3[ calculated in the CM frame. The an- 
gle of the B flight direction (0* B ) with respect to the 
beam axis provides additional discrimination since it is 



distributed as (1 — cos 2 9* B ) for B decays but flat for con- 
tinuum. The distance in the z direction (Az) between 
the signal B vertex [l3| and that of the other B is used 
in the continuum suppression if |Az| is less than 2.0 mm. 
For B events, the average value of \Az\ is approximately 
0.2 mm, whereas continuum events tend to have a com- 
mon vertex that is measured with a resolution of about 
1.0 mm. In addition, the helicity angle (0h) discrimi- 
nates between the signal and continuum events, where 
9h is the angle between the final state K + direction and 
the B meson direction in the 4> res t frame. We first cal- 
culate the individual probability density function (PDF) 
for the Fisher discriminant, co$>9* B , Az and cosOh, and 
then obtain their product, 



L 



S(qq) 



ID 



S{qq)> 



(1) 



where L l s , q -<. denotes the signal (qq) likelihood of the 
continuum suppression variable i. The PDFs for signal, 
generic B, and continuum events are obtained from the 
GEANT3-based [H| Monte Carlo (MC) simulation. The 
variable used for continuum suppression is the likelihood 
ratio (Rs) defined as 



Rs = 



Ls 



(2) 



Additional background suppression is achieved through 
the use of a _B- flavor tagging algorithm [la ], which pro- 
vides two outputs: q = ±1 indicating the flavor of the 
other B in the event, and r, which takes a value between 
and 1 and is the quality of the flavor determination. 
Events with a high value of r are considered to be well- 
tagged. The continuum background is reduced by apply- 
ing a gr-dependent selection requirement on Rs- This re- 
qurement is optimized in three qr regions for B + — ¥ 4>Tt + : 
-1 < qr < -0.5, -0.5 < qr < -0.1, and -0.1 < qr < 1. 
For B° — )• 4>tt°, since we do not distinguish the B flavor, 
we use three r intervals: < r < 0.25, 0.25 < r < 0.70, 
and 0.70 < r < 1. The requirements are chosen to maxi- 
mize a figure of merit (FOM) defined as 



FOM = 



N S 



(3) 



where Ns is the number of signal MC events in the signal 
region and N B is the number of background events esti- 
mated in the signal region by assuming B(B + — > </>7r + ) = 
2.4 x 10~ 7 and B(B° -> <jm°) = 2.8 x 10~ 7 . Our back- 
ground suppression eliminates 99.4% (99.7%) of contin- 
uum background while retaining 55.8% (43.9%) of the 
signal events for B + 07r+ (B° 4>n ). 

Backgrounds from B decays are studied using large 
MC samples. The sample size for charmless decays from 
b — > u,d, s transitions corresponds to 50 times the data 
luminosity. For B — > <f)ir + , the b — > s process B — > <pK + 
is the dominant background, arising from kaon-to-pion 
misidentification. For B° — > <fiir°, a decay with a ir° in 
the final state such as B a — > 4>Kg, is the dominant con- 
tribution. This background has a signal-like distribution 
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in Mbc- However, the background populates the nega- 
tive AE region with small overlap with the signal, so its 
contribution can be extracted from a fit. 

Signal yields for B — > 4>tt decays are obtained by per- 
forming a two-dimensional extended unbinned maximum 
likelihood (ML) fit to the observables M bc and AE. The 
likelihood is 



L 



n 



J^NiPiiM^AE^ 



(4) 



where the index i denotes signal, continuum, b — > c back- 
ground, and b — > u, d, s background components, Ni is 
the yield, Pi is the PDF for each component, and the in- 
dex j indicates the event candidate. The total signal PDF 
is described as a product of the PDFs for Mt, c and AE. 
We use the decays B + — > <pK + and B° — > D°tt° as con- 
trol samples to correct for differences between data and 
MC simulations for the fitted means and widths of Mb c 
and AE. The PDF for AE is a sum of two Gaussians for 
B + — > 4>ir + with a common mean, two widths and frac- 
tion fixed to the values obtained from a fit to B + — > <f>K + 
data, and a Crystal Ball function [TtJ with the mean and 
width hxcd to the values derived from B a — > D°ir° data 



for B Q 



The PDF for Mb c is a Gaussian func- 



tion with mean and width hxed to the values obtained 
from the respective control samples for both modes. To 
obtain the two-dimensional PDF for the continuum back- 
ground, we multiply the PDF of Afbc, for which we use 
an ARGUS [3 function, with the PDF of AE, which 
is modeled using a first-order Chebyshev polynomial for 
B+ ->■ 07T+ and B° </>vr°. Both the ARGUS shape pa- 
rameter and the AE slope are allowed to float. The PDF 
of the b — ¥ c background is modeled with two-dimensional 
histograms (2D HistoPDF) with each fixed yield derived 
from MC simulations. The b — >■ u,d, s transition back- 
grounds are modeled with two-dimensional histograms 
with fixed yields derived from MC simulations except for 
B+ -> <j)K+. The PDF for B+ -> <j)K+ is a double 
Gaussian function for AE and a Gaussian for M bc , in 
which the mean, widths, fraction and yield are fixed to 
the values derived from a fit to B + — > <fiK + data using 
the particle indentification requirement Rk,-k > 0.6 for 
kaon candidates. 

Possible backgrounds to <fi — > K + K~ decays are pre- 
dominantly from B K+K'n with / (980) ->■ K+R-, 
ao(980) — > K + K~ or a nonresonant contribution. The 
two-dimensional fit to M bc and AE alone cannot distin- 
guish the signal from other B — > K + K~tt events. We 
model B — > /o(980)7r, B — > ao(980)7r and nonresonant 
B — > K + K~it with uniform phase space distributions; 
these backgrounds are treated as additional components 
in the fits. To evaluate their contributions, we examine 
events in the </> mass sidebands, M K+K - < 1.0 GeV/c 2 
and 1.039 GeV/c 2 < M K+K - < 1.1 GeV/c 2 . We apply 
the same two-dimensional fit to the sideband events as- 
suming that signal-like events are dominated by each of 
the above three background sources. The possible con- 
tribution to the signal is then included as a background 



PDF corresponding to a signal PDF with fixed mean, 
width(s) and fraction from each component. As we can- 
not distinguish these three components, we take the non- 
resonant mode that gives the largest signal yield as the 
central value. This background contribution is found to 
be 4.7lJ'g even ts f° r B + 4 m+ an d l-6^J'g events for 
B° — > 4>tt°, derived from the data sideband. The ex- 
pected yields, 4.7 events for B + — > <f>ir + and 1.6 events 
for B° — > (f>ir°, are fixed. We summarize the PDF shape 
and expected yields (fit outputs) for various components 
in Table |U 

Figure [2] shows the AE and M^ c projections of the 
fit for the selected B candidates. There are a total 
of 373 B+ </m+ and 272 B° 0tt° candidates in 
the data sample. We determine the signal yields to 
be N S (B+ ->■ 07T+) = 4.5±|;| and N S {B° -> 07r°) = 
—2.2^2, where the quoted error is statistical only. We 
observe no significant signal for B + — > <f>Ti + or B° —> cfnr 
decays. The branching fraction B is calculated from the 
observed yield as 



B[B -> cjm) 



N B . 



><f>7T 



£data X N 



(5) 



BB 



where Ns-^^n is the signal yield, N BB - is the number 
of BB pairs (where the production rates of B + B~ and 
B°B° pairs are assumed to be equal) and edata is the 
signal reconstruction efficiency. The reconstruction effi- 
ciency is defined as 



Edata = CMC X 



^data 



£dat, 



PID 



CMC 



CMC 



PID 



(6) 



where cmc is the reconstruction efficiency from MC sim- 
ulations and the branching fractions, B(4> — > K + K~) — 
48.9% and B(n° -> 77) = 98.8%, are applied to MC 
simulations. Cdata^" {cmc Rb ) is the efficiency of the 
R s requirement from data (MC simulations), £data PID 
(eMC PID ) is the efficiency of the PID requirement from 
data (MC simulations). 

We consider the systematic uncertainties in the ef- 
ficiency, N BB and the yield extraction. The main 
sources of efficiency uncertainties are MC statistics 0.6% 
(0.8%), PID 2.0% (1.3%) and tracking 3.1% (2.0%) for 
B + -> (B° 4>n°). The uncertainty on the 

7T° efficiency is measured by comparing the yields be- 
tween T] — > 77 and 77 — > 7r 7r°7r° and is found to be 
3.0%. To evaluate the uncertainty from the efficiencies 
due to the Rs requirements, we use the control sam- 
ples B+ ->■ D°(D° ->• K+Tr-)Tt+ for B+ -> and 
B° ->■ D+(D+ -> K s n+)Tr- for B° -> ^tt . The R s un- 
certainties are 2.4% (4.1%) for B+ -> <})tt + (B° -> 0tt°). 
The uncertainty from N BB is 1.4%. The sources and 
sizes of systematic uncertainies in the signal yield ex- 
traction are listed in Table HT1 The systematic error from 
signal yield extraction is estimated by varying all fixed 
parameters by ±lcr. To obtain the errors due to the 
fixed yields of b — > u,d, s backgrounds, b — > c back- 
grounds and nonresonant B — > K + K ~tt, we vary these 
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TABLE I: Summary of the PDF's used in the measurement of B — > 4>n decays. Here CB is a Crystal Ball function and 2D 
HistoPDF is a PDF based on a histogram. Yields in the parentheses are expected values (fit outputs) for the fixed (floated) 
case. 



Mode 


B + <j>n+ 






AE M bc 


Method (Yield) 


AE M bc 


Method (Yield) 


Signal 

e + e~ — > qq process 

b — > c 

b — s> u,d,s 

B+ -> (j)K+ 

Nonresonant B —¥ K + K~ir 


Sum of two Gaussians Gaussian 
l at order poly. ARGUS 

2D HistoPDF 

2D HistoPDF 
Sum of two Gaussians Gaussian 
Sum of two Gaussians Gaussian 


Float(4.5+43) 

Float(330.0l"; 4 ) 

Fixed(7.1) 

Fixed(4.1) 

Fixed(33.8) 

Fixed(4.7) 


CB Gaussian 
1 st order poly. ARGUS 
2D HistoPDF 
2D HistoPDF 

CB Gaussian 


Float(-2.2+^) 
Float(265.6t^) 
Fixed(4.8) 
Fixed(13.5) 

Fixed (f. 6) 




FIG. 2: Projection of the data (points with error bars) in the fit region. The fit projections onto AE (left) and M^ c (right) for 
reconstructed B + — » (jm + (top) and B° — s> (jm (bottom); the sum of signal and qq (blue-dotted), qq (red-dashed), nonresonant 
B — > K + K~-k background (green-solid), other B background (magenta-solid) and the total (blue-solid). 



fixed yields by ±50%. The uncertainty from the b — > c 
backgrounds is negligible. The largest difference in yield 
between nonresonant B — > K + K~n and the other modes 
is included in the systematic error. The uncertainty from 
this difference, which is the largest contributor to the to- 
tal systematic error, is —6.3 (—2.2) events for B + — > cf)ir + 

(B '^ 07T°). 

The upper limit (Bjjl) is determined as 



C L C(B)dB 
J Q °°C(B)dB 



= 0.90, 



(7) 



where C(B) is the likelihood value and B is the branch- 
ing fraction. The branching fraction is determined as the 



number of the signal events divided by the number of 
BB pairs and the reconstruction efficiency. We include 
systematic errors by convolving the likelihood function 
with a Gaussian whose width is equal to the total sys- 
tematic error. The upper limits on the branching frac- 
tions are found to be B(B + — > 07r + ) < 3.3 x 10~ 7 and 
B(B° -> <jm°) < 1.5 x 10~ 7 at the 90% CL. The results, 
together with the central values for the branching frac- 
tions, are listed in Table ITTT1 

In summary, using 657 x 10 6 BB pairs collected at the 
T(4S) with the Belle experiment, we find no significant 
signals for B + — > </>7r + and B° —> 07r°. We set upper 
limits of B(B+ 07T+) < 3.3 xlO" 7 and 23(5° 0tt°) < 
1.5 x 10~ 7 at the 90% CL. 
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TABLE II: Summary of systematic uncertainties (events) in 
the signal yield (V) extraction. 



Sources 


B+ 07T + 


B" -»• 4>n u 


Signal PDF 


-1-0 5 

-0^6 


+0 6 

-0^4 


b — > u, d, s 


+0.1 
-0.1 


+ 0.0 
-0.1 


B + -> <t>K+ 


+ 1.8 
-1.6 




nonresonant B — ¥ K + K~ir 


±2.4 


±0.8 


Fit bias 


+0.9 
-0.0 


+0.3 
-0.0 


Peaking background modeling 


+0.0 
-6.3 


+0.0 
-2.2 


Total 


+3.1 

-6.9 


+ 1.3 

-2.4 
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TABLE III: Signal yields, measured branching fractions in- 
cluding statistical and systematic errors, and the upper limits 
including systematic uncertainties at the 90% CL. 

B + -+ (j)TT + B° -+ V 

Yield 4^1™ =2^™ 

Edata 8.4% 4.9% 

n(-\C\~ 7 \ n O+0.9+0.6 n 7 +0.6+0.4 
£S(1U ) U.8_o_ 8 _ 1 _3 -U.r_ _4_ _ 8 

B UL (W' 7 ) 3.3 1.5 
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